NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY,  CALIEORNIA 


THESIS 


MEMS-BASED  WASTE  VIBRATION  AND  ACOUSTIC 

ENERGY  HARVESTERS 

by 

Timothy  J.  Householder 

December  2014 

Thesis  Advisor: 

Dragoslav  Grbovic 

Co-Advisor: 

Bruce  Denardo 

Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


Public  repoitmg  burden  for  this  collection  of  information  is  estunated  to  average  1  hour  per  response,  mcluding  the  time  for  reviewing  instruction, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  conqileting  and  reviewing  the  collection  of  information.  Send 
comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to 
Washington  headquarters  Services,  Directorate  for  Infonnation  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 
December  2014  1 

3.  REPORT  TYPE  AND  DATES  CO^TRED 
Master’s  Tliesis 

4.  TITLE  AND  SUBTITLE 

MEMS-BASED  WASTE  VIBRATION  AND  ACOUSTIC  ENERGY 

HARVESTERS 

5.  ELTsDING  NTWIBERS 

6.  AUTHORfS)  Timothy  J.  Householder 

7.  PERFORMING  ORGANTZATION  NAME(S)  ANT)  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey.  CA  93943-5000 

8.  PERFORMING  ORGANTZATION 

REPORT  NLTMBER 

9.  SPONSORING  AIONITORING  AGENCY  NAME(S)  A.ND  .ADDRESS(ES) 
N/A 

10.  SPONSORING/MONTTORING 
AGENCY  REPORT  NT/MBER 

11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  pohcy 
or  position  of  the  Departnient  of  Defense  or  the  U.S.  Govemnient.  IRB  Protocol  number _ N/A 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  1 12b.  DISTRIBUTION  CODE 

Approved  for  public  release:  distribution  is  unlimited  |  A 

13.  AJBSTRACT  (maximum  200  words) 

Every  machine  vibrates  and  emits  noise.  This  is  unused  energy  that,  with  an  appropriate  mechanism,  can  be  reniraed 
to  the  system.  Utilizmg  an  array  of  piezoelectric  microelectromechanical  systems  (MEMS)  devices  to  harvest  tliis 
otherwise  wasted  energy,  it  is  possible  to  improve  the  efficiency  of  any  number  of  mechanical  devices. 

Piezoelectricity  is  the  mechanism  by  which  certain  crystalline  structures  generate  electric  potential  when  imder  sti'aui. 
or,  conversely,  deform  when  subjected  to  an  electric  potential.  It  is  tliis  first  effect  that  is  importarrt  to  tlris  application. 
Though  each  MEMS  device  will  generate  a  very  small  amoiuit  of  power,  a  1  m^  area  carr  contaiir  arr  array  of  millions 
of  these  devices. 

Errergy  harvesting,  corrservation.  and  efficiency  are  all  key  Department  of  Defense  (DOD)  priorities,  and  tire  imiversal 
application  of  these  devices  make  them  ideal  for  any  expeditionary  platfonn.  such  as  ships,  aircraft,  and  automobiles. 

Tliis  thesis  designs  and  tests  die  first  generations  of  acoustic  and  vibrational  piezoelectric  MEMS  devices:  including 
tinie-dependeut  finite  element  models,  microfabrication  processes,  and  the  initial  attempts  at  characterization  and 
optunization. 


14.  SL^JECT  TERMS 

MEMS.  energy  harvxsting,  piezoelectric,  resonator 

15.  NTTVIBER  OF 
PAGES 

89 

16.  PRICE  CODE 

17.  SECLRITY 

18.  SECLRITY 

19.  SECLRITY 

20.  LESHTATION  OF 

CLASSIFICATION  OF 

CLASSIFICATION  OF  THIS 

CLASSIFICATION  OF 

ABSTRACT 

REPORT 

PAGE 

ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

UU 

NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2—89) 
Prescribed  by  ANSI  Std.  239—1 8 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


MEMS-BASED  WASTE  VIBRATION  AND  ACOUSTIC  ENERGY 

HARVESTERS 


Timothy  J.  Householder 
Lieutenant  Commander,  United  States  Navy 
B.S.,  University  of  Idaho,  2002 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  2014 


Author:  Timothy  J.  Householder 


Approved  by:  Dragoslav  Grbovie 

Thesis  Advisor 


Bruee  Denardo 
Co-Advisor 


Andres  Larraza 

Chair,  Department  of  Physies 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


Every  machine  vibrates  and  emits  noise.  This  is  unused  energy  that,  with  an  appropriate 
mechanism,  can  be  returned  to  the  system.  Utilizing  an  array  of  piezoelectric 
microelectromechanical  systems  (MEMS)  devices  to  harvest  this  otherwise  wasted 
energy,  it  is  possible  to  improve  the  efficiency  of  any  number  of  mechanical  devices. 

Piezoelectricity  is  the  mechanism  by  which  certain  crystalline  structures  generate 
electric  potential  when  under  strain,  or,  conversely,  deform  when  subjected  to  an  electric 
potential.  It  is  this  first  effect  that  is  important  to  this  application.  Though  each  MEMS 
device  will  generate  a  very  small  amount  of  power,  aim  area  can  contain  an  array  of 
millions  of  these  devices. 

Energy  harvesting,  conservation,  and  efficiency  are  all  key  Department  of 
Defense  (DOD)  priorities,  and  the  universal  application  of  these  devices  makes  them 
ideal  for  any  expeditionary  platform,  such  as  ships,  aircraft,  and  automobiles. 

This  thesis  designs  and  tests  the  first  generations  of  acoustic  and  vibrational 
piezoelectric  MEMS  devices;  including  time-dependent  finite  element  models, 
microfabrication  processes,  and  the  initial  attempts  at  characterization  and  optimization. 
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I.  INTRODUCTION 


This  thesis  continues  the  work  started  by  former  Naval  Postgraduate  students 
Sarah  Gregory  and  Daniel  Hogue  in  which  the  first  generations  of  vibrational  and 
acoustic  energy-harvesting  piezoelectric  MEMS  devices  were  modeled  using  finite 
element  modeling  software,  and  the  initial  fabrication  steps  were  developed  [1].  Energy 
harvesting  devices  can  have  a  wide  range  of  designs  and  functions.  This  thesis 
investigates  piezoelectric  MEMS-based  double-ended  multifold  cantilevers. 

A,  BACKGROUND 

Developing  a  sustainable  fleet  is  at  the  forefront  of  much  of  the  Department  of  the 
Navy’s  strategic  planning.  Not  only  is  a  petroleum-dependent  fleet  expensive  and 
environmentally  burdensome,  it  also  demands  the  strategic  vulnerability  of  relying  on 
other  nations’  petroleum  exports. 

Sustainable  energy  is  among  the  most  important  issues  to  the  current 
administration  [2],  and  the  Department  of  Defense  and  its  components  have  been  directed 
“to  increase  military  effectiveness  and  mitigate  cost,  develop  more  energy-efficient 
weapons  systems,  platforms,  equipment,  and  facilities  capable  of  utilizing  multiple 
energy  sources;  invest  in  cost-effective  energy  sources  (including  alternative  energy);  and 
promote  non-materiel  and  behavior-based  solutions”  [3].  Waste  energy  harvesting  is 
identified  as  a  key  science  and  technology  objective  in  the  2012  U.S.  Marine  Corps  S&T 
Strategic  Plan  [4],  and  MEMS  energy-harvesting  devices  can  lighten  the  logistical  burden 
required  by  USMC  and  Navy  expeditionary  units. 

The  development  and  widespread  use  of  vibrational  and  acoustic  harvesting 
MEMS  would  aid  the  effort  to  meet  each  of  these  goals  in  the  Department  of  the  Navy  by 
contributing  to  the  overall  energy  portfolio  and  decreasing  the  necessity  for  use  of 
petroleum-based  generators  for  production  of  electricity  by  improving  their  overall 
efficiency. 


1 


B,  SCOPE  OF  THESIS 

Piezoelectricity  is  the  mechanism  by  which  certain  crystalline  structures  generate 
electric  potential  when  under  strain;  or,  conversely,  deform  when  subjected  to  an  electric 
potential.  Simulations  performed  by  Gregory  and  Hogue  [1]  show  it  is  possible  to  harvest 
waste  energy  to  power  small  electronic  devices  by  developing  a  piezoelectric 
microelectromechanical  system. 

This  thesis  is  primarily  focused  on  the  completion  of  the  microfabrication  steps 
initially  developed  by  Gregory,  characterization  of  the  completed  devices,  and  feedback 
of  the  obtained  parameters  (previously  unavailable)  into  the  finite-element  modeling 
simulations  conducted  by  Hogue  in  order  to  obtain  a  more  accurate  model.  Other  avenues 
of  investigation  include  reactive  ion  sputtering  of  aluminum  nitride  (AIN)  using  an 
aluminum  target  with  nitrogen  gas,  the  use  of  SU-8  photoresist,  vice  silicon  (Si),  as  a 
structural  layer,  and  analyzing  the  effect  of  additional  multifold  legs  on  resonance 
frequencies,  all  in  an  effort  to  bring  a  complete  microfabrication  process  to  NPS. 

The  remaining  microfabrication  process  steps  have  been  completed  to  finish  the 
device,  as  well  as  some  process  refinement.  Experimental  results  have  been  fed  back  into 
the  model  developed  in  [1]  in  order  to  refine  and  optimize  the  model.  These  MEMS 
devices  may  be  tuned  for  60Hz  or  50Hz  to  return  waste  energy  from  El.S.  El.S.  or  foreign 
electrical  distribution  systems,  or  can  be  designed  for  a  swath  of  resonance  frequencies  in 
order  to  design  the  devices  for  specific  vibration  instruments,  or  for  bulk  energy 
harvesting. 
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II.  PHYSICS  WITHIN  THE  DESIGN 


Energy-harvesting  devices  can  have  a  wide  range  of  designs  and  functions.  This 
chapter  introduces  MEMS  and  describes  basic  principles  of  piezoelectric  energy 
harvesting  devices. 

A,  MICROELECTROMECHANICAL  SYSTEMS 

hollowing  closely  behind  the  semiconductor-based  integrated  circuit  industry,  the 
MEMS  field  has  capitalized  upon  the  significant  advances  made  in  photolithography  and 
Si-based  micro  fabrication  [5].  MEMS  are  in  use  today  in  many  systems,  such  as  motion- 
sensors  in  mobile  phones  and  gaming  devices,  pressure  sensors  in  automobile  tires  and 
airbags,  and  a  multitude  of  small-scale  devices  such  as  microphones  and  infrared 
cameras. 

Microelectromechanical  systems,  due  to  their  small  size  and  low  cost,  are  ideal 
for  bulk  energy-harvesting  applications.  The  dimensions  of  a  typical  MEMS  device  are 
on  the  order  of  microns,  so  one  array  can  consist  of  millions,  perhaps  billions  of  MEMS 
devices.  While  each  individual  device  may  not  provide  much  usable  energy,  an  array  of 
MEMS  devices  can  have  significant  returns  with  a  very  small  footprint. 

B,  PIEZOELECTRIC  EFFECT 

Piezoelectricity  is  highly  scalable  and  can  provide  a  very  efficient  means  of 
energy  conversion  [6].  Because  of  this,  piezoelectric  devices  have  become  prevalent  in 
modem  equipment,  from  flame  lighters  to  MEMS  devices  to  ultrasonic  transceivers  for 
material  testing  [7].  Eor  a  material  to  be  piezoelectric  it  must  be  a  noncentrosymmetric 
crystal;  that  is,  it  must  have  no  center  of  symmetry. 

1.  Theory 

When  a  stress  (force  per  unit  cross-sectional  area)  is  applied  to  a 
noncentrosymmetric  crystal  in  certain  directions,  it  causes  a  strain  (fractional  change  in 
length),  which  shifts  the  center  of  mass  of  the  respective  opposing  ions.  This  shift  can 
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create  an  electric  dipole,  which  is  not  necessarily  along  the  same  direction  as  the  applied 


stress. 


If  the  piezoelectric  coefficient  is  expressed  as  a  matrix  there  is  a  linear 
relationship  between  the  applied  mechanical  stress,  7},  and  the  induced  polarization.  Pi,  as 
well  as  a  linear  relationship  between  the  applied  electric  field  Et,  and  the  indueed  strain, 
Sj.  The  indices  i  and  j  indicate  direction,  known  as  the  ‘hy-mode.”  Many  piezoelectric 
materials  have  more  than  one  possible  resonanee  mode,  describing  the  differenee 
between  the  directions  of  applied  stress  and  polarization,  which  is  determined  by  the 
crystalline  geometry  [7].  Note  the  piezoeleetric  coefficient  dy  is  the  same  in  both 
equations  below,  which  is  due  to  reciprocity: 


=  dyE, 


(1) 


Piezoelectrics,  like  all  erystals,  are  frequently  identified  by  their  axes:  a  and  b 
(horizontal)  or  c  (vertical).  Figure  1  shows  stresses  on  the  legs  of  a  piezoelectrie  device. 


Eigenfrequency=  366, 472229+2. 6760471 
Surface:  Stress  tensor,  xy  conoponent  (N/m^) 


Figure  1.  Computer-generated  finite-element  model  of  a  piezoeleetrie 
device,  showing  stress  magnitude  by  color  during  deformation. 


2,  Energy  Conversion  and  Resonance  Frequency 

The  effectiveness  of  the  energy  eonversion  process  (not  to  be  confused  with  the 
efficiency)  can  be  deseribed  with  an  electromechanieal  eoupling  factor 
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^2  _  Energy  converted 
Energy  input 
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The  piezoelectric  portion  of  any  circuit  will  oscillate  and,  due  to  its 
electromechanical  coupling,  can  be  modeled  as  an  LRC  “tank”  circuit  with  the  equivalent 
formula  for  resonance  frequency.  As  shown  in  [6],  a  piezoelectric  energy  conversion 
device  has  an  efficiency  of: 


2 


2 


r  e  ] 


(3) 


2 

where  k  is  the  electromechanical  coupling  factor  and  Q  is  the  quality  factor.  The  basic 
assumptions  used  in  this  derivation  are  that  the  maximum  efficiency  will  occur  at 
resonance  frequency  [6],  and  that  the  device  can  be  modeled  as  the  equivalent  circuit 
shown  in  Figure  2. 


Figure  2.  LRC  circuit  model  of  a  piezoelectric  device. 


where  L  represents  mass,  C  is  stiffness,  R  is  damping,  and  Zl  is  the  load  in  parallel  with 
the  piezoelectric  device’s  electrodes.  The  resonance  frequency  is 
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An  unavoidable  side  effeet  of  the  physieal  use  of  sueh  a  deviee  is  the  capaeitance 
of  the  eleetrodes  in  parallel  with  the  erystal,  Co.  This  causes  an  “antiresonance 
frequency.” 
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Just  as  with  any  other  LRC  circuit,  piezoelectric  crystals  have  a  quality  factor  Q 
that  can  be  described  in  the  same  way  [6]; 


iTTf^RC  ■ 


(V) 


This  analysis  becomes  ungainly  when  dealing  with  a  more  complicated  design,  or 
when  analyzing  an  array  of  piezoelectric  devices.  This  thesis  uses  finite-element- 
modeling  software  for  in-depth  analysis. 


C.  GENERATING  ELECTRICAL  CURRENT 

Power  is  the  rate  of  energy  transfer,  and  in  electrical  systems  is  the  product  of 
current  and  voltage.  In  order  for  an  energy-harvesting  device  to  be  useful  it  must  have 
sufficient  current  capacity  to  deliver  sufficient  power  to  its  load. 

Piezoelectric  devices  can  generate  high  electric  potential  differences  (voltage); 
however,  the  current  capacity  of  each  device  is  severely  limited.  Finite  element  models 
from  [1]  suggest  that  an  ideal  load  exists  to  which  a  single  device  can  deliver  the 
maximum  power,  as  seen  in  Figure  3.  This  analysis  will  be  addressed  in  Chapter  VIII. 
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Figure  3.  Time-dependent  eleetrieal  power  output  of  device  across  varying 
load,  driven  at  0.1  mm  displacement,  from  [1]. 

Moreover,  MEMS  devices  can  be  connected  in  series  to  step-up  the  voltage  and  in 
parallel  to  step-up  the  current  capacity  of  the  power-generating  array.  This  analysis  will 
be  addressed  in  Chapter  VIII. 
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III.  MATERIALS  SELECTION  AND  ANALYSIS 


The  careful  selection  of  materials  is  necessary  to  ensure  an  operational,  reliable, 
and  cost-effective  device.  This  chapter  describes  this  selection  process  and  evaluates 
lower-cost  alternatives. 

A,  ALUMINUM  NITRIDE 

Aluminum  nitride  is  a  common,  inexpensive  piezoelectric  material  with  a  high 
piezoelectric  modulus;  that  is,  sufficient  potential  is  generated  for  a  given  strain.  It  is  well 
suited  to  thin-film  deposition  and  micromachining  through  wet  etching,  and  is  thus  useful 
for  MEMS  applications. 

The  preferred  modes  of  electromechanical  coupling  for  AIN  are  the  33  (parallel) 
and  31  modes  (perpendicular).  For  the  purposes  of  this  thesis  a  highly  c-axis  oriented 
AIN  layer  that  excites  the  31  mode  is  preferred  [1].  There  are  two  basic  ways  to 
determine  orientation;  examining  morphology  with  a  scanning  electron  microscope 
(SEM),  and  analyzing  X-ray  (Bragg)  diffraction  (XRD),  as  shown  in  Figure  4. 


3272  nm 


4/16/2013  HV  mag  o  WD 
630:03AM  2  00kV  SOOOOOx  172m 


Position  Y“10\ 

Figure  4.  a)  SEM  and  b)  XRD  of  c-axis  oriented  thin-film  AIN,  from  [1]. 
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B,  BIMATERIAL 

Finite  Element  Modeling  shows  that  a  thin-film  piezoelectric  layer  can  generate  a 
greater  potential  difference;  however,  this  can  lead  to  a  structurally  unsound  device.  A 
structure  layer  is  necessary  in  order  to  have  a  device  that  can  oscillate  without  failure. 
Silicon,  due  to  its  relative  abundance  and  affordability,  is  among  the  most  common 
materials  used  for  this  purpose. 

C,  SUBSTRATE 

Due  to  their  relative  abundance  and  common  use  in  CMOS  circuits,  we  use  Si 
substrates.  The  most  effective,  and  most  expensive,  is  Silicon-on-insulator  (SOI).  Other 
designs  use  deposited  layers  atop  a  bare  silicon  wafer. 

SOI  wafers  typically  have  a  400  pm-600  pm  Si  substrate  with  a  thin  layer  (0.5- 
10  pm)  of  silicon  oxide  (SiOi)  between  the  substrate  and  the  surface  layer  of  1  pm- 
20  pm  Si.  This  layout  easily  lends  itself  to  the  use  of  Si  as  a  structure  layer. 

Other  materials  are  being  investigated  for  use  as  a  structure  layer,  including 
aluminum  and  SU-8  Photoresist.  The  completion  of  this  investigation  is  left  for  future 
work  (see  Section  VIII. B). 

D,  DEVICE  ANALYSIS 

The  finite-element-modeling  software  COMSOL  is  used  to  perform  computer 
analysis  of  our  design  prior  to  fabrication.  The  physical  device  is  tested  and  the 
experimental  results  are  fed  back  into  the  computer  model.  This  revises  the  given 
parameters  to  modify  and  update  computer  models  to  more  accurately  reflect  reality. 
After  sufficient  accuracy  of  the  model  is  achieved,  multiple  iterations  of  new  generations 
of  the  energy  harvesting  devices  can  be  optimized  and  designed  without  the  need  for 
time-consuming  and  costly  microfabrication  and  characterization.  This  will  lead  to  the 
ability  to  microfabricate  highly  optimized  devices.  This  iterative  process  is  left  for  future 
work  (see  Chapter  VIII). 
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IV.  DESIGN 


There  are  several  competing  priorities  for  design  criteria.  The  structure  must  be 
robust  enough  to  withstand  cyclic  loading  with  pressures  and  displacements  seen  in  real- 
world  applications.  A  U.S.  Navy  gas  turbine  generator  can  vibrate  up  to  6  thousandths  of 
an  inch,  or  150  pm,  before  alarming  [8].  At  its  ideal  operating  speed  of  3600  RPM,  a 
vibrational  energy  harvesting  device  must  withstand  several  billions  of  cycles  each  year 
without  failure. 

The  device  must  have  a  resonance  frequency  that  is  low  enough  for  real-world 
applications.  For  use  with  U.S.  electrical  generators  this  frequency  is  60  Hz,  for 
multinational  generators  it  is  50  Hz  (see  paragraph  II.B.2). 

The  device  must  also  have  good  electrical  properties:  it  must  induce  a  high 
electric  potential,  and  have  a  high  current  capacity  (see  sections  II.B  and  II. C). 


A.  SUSPENDED,  MULTIFOLD  MICROBRIDGE  DESIGN 

MEMS  devices  tend  to  have  resonance  frequencies  in  the  kHz  and  MHz  ranges, 
though  most  useful  applications  are  in  the  Hz  ranges.  The  fundamental  mode  resonance 
frequency  for  a  MEMS  double-ended  cantilever  device  with  a  central  point  load  F  can  be 
approximated  with 


f 


13.86  \Elg 


In 


(8) 


and 


1  = 


(9) 


where  /  is  the  moment  of  inertia,  E  is  Young’s  modulus,  g  is  the  acceleration  of  gravity 
near  the  surface  of  the  earth,  and  I,  w,  &  t  are  the  device’s  length,  width,  and  thickness, 
respectively  [5].  In  order  to  design  a  device  with  a  resonance  frequency  in  usable  ranges, 
the  device  must  be  long  and  thin,  with  a  majority  of  its  mass  in  the  center.  The  most 
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compact  possible  geometry  is  used:  long,  multifold  legs  (providing  large  effeetive  length) 
with  a  massive  eenter  pad. 


B,  GENERAL  DIMENSIONS  AND  THEIR  IMPLICATIONS 

The  highest  eleetrie  potential  will  be  indueed  at  the  point  of  the  greatest  stress,  as 
shown  in  Equation  (1).  The  point  of  greatest  stress  will  also  be  the  most  suseeptible  to 
failure. 


The  deviee  legs  must  be  thiek  enough  and  wide  enough  to  eompensate  for  the 
longer  legs  and  withstand  the  repetitive  eyelie  stresses  of  vibration  without  suffering 
failure  or  early  fatigue.  For  a  double-ended  eantilever  (fixed-fixed  beam)  the  maximum 
stress  will  oeeur  in  the  middle  of  eaeh  leg,  whieh  ean  be  modeled  as  a  fixed-free  beam 
under  point  loading,  as  follows: 


e 


2EI  ’ 


(10) 


where  6  is  the  defieetion  angle,  whieh  eorrelates  to  stress  [5].  Note  the  similarity  to 
Equation  (8).  A  low  /  or  high  t  is  desired  to  reduee  the  maximum  stress  felt  by  the 
eantilever  strueture.  This  is  in  direet  eonfliet  with  the  requirements  for  a  low  resonanee 
frequeney  and  a  high  indueed  eleetrieal  potential. 


C.  SPECIFIC  DESIGNS 

This  thesis  eovers  multiple  lines  of  investigation,  eaeh  using  a  different  design. 
The  design  used  in  [1]  is  awaiting  Deep  Reaetive  Ion  Etehing,  its  final  fabrieation  step, 
prior  to  experimental  analysis.  A  set  of  similar  designs  was  used  to  develop  an  ad-hoe 
formula  for  resonanee  frequeney  vs.  eantilever  leg  length,  speeifieally  with  regard  to  the 
number  of  folds  in  the  multifold  legs;  these  designs  are  also  intended  to  be  used  in  the 
examination  of  A1  and  of  SU-8  as  potential  strueture  layers.  The  eontinued  analysis  of 
these  designs  is  left  for  future  work  (see  Chapter  VIII). 

Finally,  a  more  robust  design  was  developed  for  fabrieation  by  a  third-party 
eommereial  foundry,  MEMSCap  [1,  9].  A  new  method  of  fabrieating  piezoeleetrie 
MEMS  deviees,  known  as  PiezoMUMPs,  is  under  development  and  this  thesis  provides 
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some  data  to  help  characterize  the  PiezoMUMPs  process  [1].  The  experimental  data  in 
Chapter  VII  come  from  the  PiezoMUMPs  devices. 

Each  design  was  modeled  with  COMSOL  and  designed  using  MEMS  design 
software  called  MEMSPro. 

1,  Design  1:  PiezoMUMPs  Fabrication 

The  first  generation  PiezoMUMPs  design  in  [I]  did  not  survive  shipping  and  thus 
could  not  be  tested.  The  process  itself  was  not  the  problem:  the  final  step  of  cleaving  the 
wafer  into  separate  die  while  on  sticky  tape  prior  to  shipping  was  where  this  issue 
occurred.  The  maximum  displacement  of  the  center  pad  would  exceed  the  substrate 
thickness  and  stick  to  the  shipping  tape.  During  the  attempted  removal  of  the  device  from 
the  shipping  tape,  the  center  pad  (along  with  most  of  the  legs)  would  remain  affixed  and 
shear  off.  Additionally,  both  electrode  tabs  tapped  off  the  top  of  the  legs  and  thus  would 
not  have  obtained  the  optimum  electric  potential  generated  from  the  top  to  the  bottom  of 
the  device. 

The  second  generation  device  was  designed  with  thicker  legs  to  reduce  the 
maximum  displacement  of  the  center  pad,  and  with  electrode  pads  placed  to  obtain  the 
optimal  “top-to-bottom”  electric  potential.  The  tradeoff  with  this  design  is  a  higher 
theoretical  resonance  frequency.  The  center  pad  continued  to  stick  to  the  shipping  tape. 

The  third  generation  device  is  more  robust  device  with  fewer  legs;  connected  at 
the  midpoint,  rather  than  the  corner,  of  the  center  pad;  and  with  temporary  “bridges”  at 
the  corners  that  pin  the  center  pad  to  the  substrate  during  shipping  (Eigure  5).  This  device 
has  a  higher  theoretical  resonance  frequency  than  the  two  previous  generations. 

These  three  PiezoMUMPs  generations  are  shown  in  Eigure  5. 
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Figure  5. 


PiezoMUMPs  generations.  Note  the  bridges  attaehed  to  the  eomers 
of  the  eenter  pad  in  the  third  generation  deviee. 


2,  Design  2:  NPS  Fabrication  with  SOI 

There  have  been  no  changes  to  the  designs  for  NPS  fabrication  used  in  [1].  This 
design  will  test  2-,  4-,  and  6-legged  structures.  Figure  6  shows  the  dimensions  of  the  6- 
legged  structure  on  the  photolithography  mask,  generated  with  MEMSPro,  of  this  device. 
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Figure  6.  Graphic  of  the  mask  for  a  6-legged  piezoelectric  energy  harvester, 
annotated  with  structure  notes  and  dimensions,  from  [1],  The 
device  is  grey,  the  substrate  is  light  blue,  and  the  trench  is  white. 


This  design  uses  SOI  wafers  which,  though  more  costly,  enable  full  production  at 
NFS.  Since  the  pre-fabricated  SOI  wafer  already  contains  the  structure  layer,  only  two 
photolithography  masks  are  required  for  complete  microfabrication,  the  piezo/structure 
mask,  and  the  backside  (trench)  mask.  This  process  will  be  described  in  Chapter  V.  Due 
to  the  higher  cost  of  SOI  wafers  (5-7  times  the  price  of  a  bare  silicon  wafer),  however, 
alternative  structure  materials,  such  as  SU-8  and  aluminum,  are  also  being  investigated. 

3,  Design  3:  NFS  Fabrication  with  SU-8 

If  structurally  viable,  SU-8  has  some  beneficial  characteristics:  it  is  quite 

inexpensive  and  it  can  be  deposited  quickly  and  easily  via  spin-coating  in  any  laboratory. 

Due  to  its  low  Young’s  modulus  (as  low  as  2  GPa  [10]  and  as  high  as  5  GPa  [11])  lower 

resonance  frequencies  and  greater  stress  to  the  piezo  layer  can  be  achieved  at  much 
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smaller  device  size  than  would  be  possible  with  silicon-based  devices,  allowing  for  more 
dense  packing  and  thus  a  larger  number  of  devices  per  unit  area. 

Two  separate  wafer  designs  have  been  prepared  to  test  SU-8  as  a  structure  layer. 
Design  2:  NFS  Fabrication  with  SOI  is  used  for  a  known  comparison,  and  the  second 
design  will  test  the  frequency  response  of  SU-8  over  a  larger  range  of  leg  length, 
identified  by  the  number  of  folds  per  side,  as  shown  in  Figure  7. 


Figure  7.  MEMSPro  images  of  SU-8  designs,  1-  to  8-legged  structures. 


Since  the  piezo  layer  is  deposited  directly  atop  the  structure  layer,  only  two 
photolithographic  masks  are  required  for  complete  microfabrication,  just  as  with  Design 
2:  NPS  Fabrication  with  SOI.  This  process  will  also  be  described  in  Chapter  V.  These 
processes  can  be  greatly  simplified  by  expanding  AIN  reactive  sputtering  capabilities  at 
NPS. 
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V.  MICROFABRICATION  PROCESS 


As  explained  in  Chapter  IV,  there  were  multiple  mierofabrieation  teehniques 
evaluated.  This  ehapter  deseribes  refinement  of  the  NFS  mierofabrieation  proeess. 

A,  GENERAL 

1,  Sputter  Deposition  of  AIN 

Sputter  deposition  of  AIN  using  faeilities  at  NFS  proved  diffieult.  The  instrument 
is  intended  to  eoat  thin  films  of  just  a  few  nanometers  for  eleetron  mieroseopy,  and  this 
eapability  did  not  seale  to  deposition  thieknesses  of  several  hundreds  of  nanometers  onto 
a  4-inoh  wafer.  The  instrument  was  also  not  designed  for  use  with  more  than  one  gas, 
intended  to  be  argon.  Finally,  the  sputter  eoater  at  NFS  did  not  have  any  temperature 
eontrol. 

Initial  attempts  to  sputter  AIN  onto  a  Si  substrate  examined  varying  ratios  of 
argon  (Ar)  and  nitrogen  (N2)  gas,  using  a  Cressington  Sputter  Coater  208HR.  Literature 
shows  better  e-axis  orientation  using  higher  temperature,  lower  pressure,  more  power, 
and  greater  pereentage  of  N2  in  the  Ar/N2  blend.  The  only  tradeoff  was  in  the  Ar/N2 
blend;  more  Ar  meant  faster  sputtering  and  less  residual  stress  [12-19].  Different  Ar/N2 
blends  were  tested  on  eaeh  of  four  samples  taken  from  the  same  Si  wafer.  A  99.99%  pure 
A1  target  was  used  for  all  reaetive  sputtering.  Table  1  summarizes  the  experimental 
parameters,  and  Figure  8  shows  the  atomie  speetroseopy  results. 
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Table  1. 


Summary  of  attempted  reaetive  sputtering  lOnm  AIN,  with  varying 

Ar/N2  blend. 


Sample 

# 

N2/(N2-Ar) 

[%] 

Chamber  Pressure 

[mBarr] 

Sputter  Current 

[mA] 

1 

Control 

(Not 

Sputtered) 

2 

0 

0.02  ±0.01 

78  ±2 

3 

100 

0.03  ±0.01 

76  ±2 

4 

50 

0.02  ±0.01 

76  ±3 

Figure  8.  Atomic  spectroscopy  comparison  of  sputtering 
with  varying  Ar/N2  blends. 


The  plasma  color  indicates  the  amount  of  Ar  and  N2,  as  seen  in  Figure  9.  The 
presence  of  oxygen  also  appeared  to  provide  a  more  reddish  hue. 
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Figure  9.  Sputtering  of  A1  target  to  Si  substrate  with  a)  100%,  b)  50%,  and  c) 

0%  N2. 


Figiue  10  shows  atomic  spectroscopy  of  an  attempted  deposition  of  AIN  using 
100%  Ni  with  a  fiill  Si  wafer  substrate.  The  high  concentration  of  oxygen  indicates  that 
ahuuina  (AI2O3)  was  preferentially  deposited,  followed  by  a  roughly  equivalent  layer  of 
AIN. 


Figiue  10.  Atomic  spectroscopy  of  sputtered  AI2O3/AIN  using  100%  N2. 

Depeudmg  on  crystalline  orientation,  the  Si  substrate  has  a  native  oxide  layer 
approximately  1  nm  thick  [20],  and  the  A1  target  has  a  native  oxide  layer  approximately 
5  nm  thick  [21].  Table  2  shows  roughly  equal  concentrations  by  number  of  oxygen  and 
nitrogen,  indicating  gr  owth  of  both  AI2O3  and  AIN. 
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Table  2.  Atomic  composition  of  sputtered  AI2O3/AIN. 


Element 

Wt% 

At% 

N 

08.68 

14.88 

0 

09.65 

14.48 

A1 

23.34 

20.77 

Si 

58.33 

49.87 

The  morphological  properties  in  Figure  1 1  appear  very  similar  to  AIN  and  show 
the  desired  “pebble  type”  c-axis  orientation.  These  data  are  consistent  with  an  initial 
deposition  layer  of  AI2O3  on  top  of  the  natural  Si02  layer;  then,  when  sufficient  oxygen 
inventory  is  reduced,  AIN  is  deposited  on  top. 


Filmetrics  measurements  show  a  roughly  Gaussian  distribution  of  sputtered  layer 
thickness  atop  a  1 0-cm  wafer,  varying  from  approximately  90  nm  at  the  edges  to  200  nm 
in  the  center.  The  concentric  rings  in  Figure  12  indicate  the  non-uniform  deposition 
thickness. 
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Figure  12.  Whole  wafer  after  attempted  AIN  deposition  at 

Naval  Postgraduate  School. 


B.  INDIVIDUAL  MICROFABRICATION  SEQUENCE 

This  thesis  further  explores  the  NPS  fabrication.  These  steps  were  developed  in 
[1]  and  are  listed  again  here,  with  some  clarification  and  updates,  for  continuity: 

Fabrication  begins  with  an  SOI  wafer,  a  section  of  which  is  shown  in 
Figure  13a.  Then,  fabrication  materials  are  deposited  and  etched  in  the 
following  order: 

1.  AIN  is  deposited  by  AC  powered  S-gun  magnetron  for  AIN  reactive 
sputtering,  as  shown  in  Figure  13b. 

2.  The  wafer  is  primed  by  spinning  the  wafer  at  3000  RPM  under  several 
milliliters  of  MCC  Primer,  an  organic  solvent  solution,  to  prepare  the 
wafer  for  photoresist  application. 

3.  7  microns  of  SPR  220  photoresist  are  spin-coated  on  the  primed  wafer 
at  3000  RPM  for  45  s.  The  wafer  is  then  soft  baked  at  1 15°C  for  90  s. 

4.  Next,  the  photoresist  is  exposed  under  the  mask  described  in  Chapter 
IV. C. 2  at  between  240  and  250  Integra,  with  a  print  gap  of  30  microns. 

5.  The  wafer  rests  for  30  minutes  to  allow  for  more  solvent  in  the 
photoresist  to  evaporate.  After  at  least  30  minutes,  the  wafer  is  post¬ 
exposure  baked  at  1 15°C  for  90  s. 

6.  Photoresist  is  developed  for  1  minute  in  CD-26  developer,  to  remove 
exposed  photoresist. 

7.  The  AIN  is  isotropically  etched,  using  85%  H3P04  at  85°C  for  60-74  s, 
as  shown  in  Figure  13c. 

8.  Using  the  AIN  as  a  mask,  a  reactive  ion  etch  will  remove  the  top  layer 
of  Si,  defining  the  structure,  as  shown  in  Figure  13d. 
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9.  The  top  AIN  and  Si  structures  are  protected  by  spinning  the  wafer 
surface  with  SPR  220-7,  and  then  baking  for  four  hours  at  90°C  [20]. 

10.  The  bottom  side  of  the  wafer  is  primed,  as  it  was  in  step  2.  Next,  it  is 
spin  coated  with  SPR  220-7  under  the  same  conditions  in  step  3. 

11.  The  trench  mask  is  aligned,  using  a  backside  aligner,  which  is 
available  at  Stanford  Nanofabrication  Facility  at  Stanford  University. 

12.  Expose  the  wafer  at  220  Integra,  with  a  print  gap  of  30  microns. 
Develop  for  3  minutes  in  CD-20.  The  result  is  shown  in  Figure  13e. 

13.  Etch  the  trench  through  the  bottom  Si  layer  using  deep  reactive  ion 
etcher. 

14.  Remove  the  photoresist  from  the  bottom  layer,  by  soaking  in  acetone 
for  5  minutes  and  rinsing  with  isopropanol. 

15.  Use  reactive  ion  etch  on  the  Si02  to  complete  the  trench  and  release 
the  energy  harvesting  structure,  as  shown  in  Figure  13f 


Figure  13.  The  Naval  Postgraduate  School  microfabrication.  Silicon  is  shown 
in  grey.  SiOi  is  shown  in  green.  AIN  is  shown  in  blue,  from  [1]. 
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The  following  additional  steps  to  ereate  an  SU-8  strueture  layer  oeeur  prior  to  step 
1,  above. 

1 .  Starting  with  a  bare  Si  wafer,  500  nm  of  SU-8  is  spin-eoated  at  3000  RPM 
and  then  softbaked  for  60  seeonds  at  95°C. 

2.  The  SU-8  is  exposed  without  a  mask  at  595  Integra,  then  post-exposure 
baked  for  120  seconds  at  95“C,  and  finally  hardbaked  for  30  minutes  at 
250°C. 

The  wafer  is  then  treated  like  an  SOI  wafer  for  the  purposes  of  the  above 
procedure  without  the  need  for  step  15. 
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VI.  MATERIAL  AND  DEVICE  CHARACTERIZATION 


Characterization  of  each  wafer  and  each  device  oeeurred  with  the  Scanning 
Electron  Microscope,  and  with  the  various  devices  in  the  Clean  Room. 

A,  ELECTRON  MICROSCOPY 

SEM  images  were  taken  at  various  loeations  on  the  seeond  generation 
PiezoMUMPs  deviee.  The  first  image  in  Eigure  14  shows  the  entire  deviee  with  the 
eenter  pad  and  three  legs  stuek  to  the  paeking  tape.  Eaeh  image  shows  lateral  separation 
between  the  Si  structure  layer  and  the  AIN  piezo  layer.  Images  of  the  center  pad  also 
show  lateral  separation  between  the  AIN  piezo  layer  and  the  A1  pad  layer.  The  even 
horizontal  striations  on  the  structure  layer  are  indicative  of  the  deep  reactive  ion  etching 
(also  called  BOSCH)  proeess  and  the  slight  irregularities  on  the  piezo  and  pad  layers  are 
indicative  of  wet  etehing  [22].  Note  the  e-axis  orientation  of  the  AIN  shown  at  the  top  of 
the  seeond  image  from  a). 
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Figure  14.  SEM  images  of  seeond  generation  PiezoMUMPs  deviee. 


B,  PHYSICAL  AND  OPTICAL  MEASUREMENTS 

Thiekness  measurements  were  eondueted  using  the  following  deviees  in  the  Clean 
Room;  KLA  D120  Contaet  Profilometer,  a  Pro-4  4-point  eontaet  probe,  Filmetrics 
speetroseopic  thiekness  measurement  instrument,  and  Zygo  NanoView  7100  Optieal 
Profilometer. 
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The  optical  profilometer  was  used  to  measure  the  NPS  Fabricated  on  SOI  device, 
shown  in  Figure  15.  The  height  of  the  device  is  uniform  across  the  center  pad  and  its  legs, 
approximately  2.65  pm. 
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Figure  16.  Optical  profilometer  measurements  of  device.  The  curvature  of  the 

center  pad  and  legs  is  expected  and  is  due  to  residual  stress  from 
the  microfabrication  process.  (Note:  scale  is  exaggerated  in  z 

direction). 
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C.  LAYER  THICKNESS 

The  best  values  for  use  in  finite-element  analysis  were  determined  above  and  are 
listed  in  Table  3. 


Tables.  Deviee  thiekness. 


Structure  (Si) 

Piezo  (AIN) 

Pad  Metal  (Al) 

Nominal  (Design) 

10.0  pm 

500  nm 

1 .02  pm 

Optical  Profilometer 

10.1  pm 

498.9  nm 

0.99  pm 

Filmetrics 

- 

493.2  nm 

- 

Best  Value 

10  pm 

493  nm 

1  pm 

D,  NANOINDENTER  (YOUNG’S  MODULUS) 

A  die  with  a  broken  device  from  the  third  generation  PiezoMUMPs  device  was 
analyzed  using  an  Agilent  Nanoindenter  G200  to  find  Young’s  modulus  of  the  Si  and  the 
AIN.  Young’s  modulus,  or  the  modulus  of  elasticity,  is  a  proportionality  constant  that 
relates  an  applied  stress  to  the  strain  a  given  material  will  display  in  its  elastic  region. 
This  can  be  compared  to  equations  (1),  (8),  and  (10),  and  is  often  expressed  as 

(j  =  Es,  (11) 

where  E  is  Young’s  modulus  (as  in  equations  (8)  and  (10),  not  to  be  mistaken  for  the 
induced  electric  field,  Ei  in  Equation  (1)),  cr  is  the  stress  {Sj  in  (1)),  and  5  is  the  strain  (7) 
in  (1))  [5].  Since  Young’s  modulus  is  the  measure  of  a  material’s  elasticity,  it  must  be  as 
accurate  as  possible  in  order  to  develop  computer  models  that  can  accurately  determine 
resonance  frequency,  displacement,  and  stresses  on  a  device. 

The  Nanoindenter  impresses  a  triangular  diamond  tip  into  a  material  (see  Figure 
17  )  to  determine  the  Young’s  modulus  of  bulk  material,  given  Poisson’s  ratio  (v);  and  of 
thin  films,  given  the  film  thickness,  Poisson’s  ratio,  and  Young’s  modulus  of  the 
substrate.  The  thin  film  analysis  mathematically  removes  the  expected  effect  of  the 
substrate  on  the  thin  film  results.  This  thesis  uses  the  following  values  for  Poisson’s  ratio: 
vsi  =  0.28  [23],  and  vain  =  0.22  [24];  and  anticipates  the  following  values  for  Young’s 
modulus:  Esi  =180  GPa  [23],  and  Fain  =  331  GPa  [24]. 
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Figure  17.  Ten  nano  indentations  of  AIN  film. 


Using  vsi  =  0.28,  ten  500  nm  deep  test  points  were  taken  on  the  silicon  layer;  one 
(test  #4)  was  anomalous.  Results  are  shown  in  Figure  18.  After  approximately  50  nm  the 
value  is  consistent  and  has  a  very  low  slope.  Averaging  the  Young’s  modulus  values 
between  100  nm  and  200  nm,  the  calculated  value  is  ^si  =  173.3  GPa,  with  standard  error 


SE, 


Sample 


^  ^ m  "^Sample  ) 


[(n-l) 


(12) 


where  5  is  the  sample  standard  deviation,  is  the  value  of  an  individual  test,  £'sampie  is 
the  sample’s  average  value,  and  n  is  the  number  of  tests.  For  this  run  of  tests,  SE^i  =  0.8 
GPa  (roughly  5%). 
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silicon 


Figure  18.  Nanoindenter  results  for  Young’s  modulus  of  silicon,  9  of  10  tests. 

Esi  =173.3  GPa,  =  0.8  GPa. 


Using  Vain  =  0.22,  fdm  thickness  of  493  nm,  and  substrate  (Si)  data  from  above, 
twenty-four  400  nm  deep  test  points  were  taken  on  the  aluminum  nitride  thin  film  layer; 
two  (test  #’s  11  &  16)  were  anomalous.  Results  are  shown  in  Figure  19.  Above  150  nm 
the  slope  becomes  fiat;  however,  the  variance  between  individual  tests  grows 
substantially  after  300  nm.  Averaging  the  Young’s  modulus  values  between  35%  and 
60%  of  film  thickness  (173  nm  to  296  nm),  the  calculated  value  is  ^ain  =  384.0  GPa,  with 
standard  error  3.8  GPa  (roughly  10%). 
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Aluminum  Nitride 


Figure  19.  Nanoindenter  results  for  Young’s  modulus  of  aluminum  nitride,  22 

of  24  tests,  ^ain  =  384.0  GPa,  SEaw  =  3.8  GPa. 


It  should  be  noted  that  ^ain  has  excessive  curvature  between  50  nm  and  100  nm. 
This  is  not  observed  in  the  ^si  data  in  Figure  18,  or  in  the  fused  silica  reference  standard. 
Though  there  is  greater  variance  between  50  nm  and  100  nm,  the  reference  results  are 
consistent.  This  implies  the  low-depth  values  for  AIN  represent  some  kind  of  surface 
affect.  Further  investigation  of  this  phenomenon  may  improve  computer  models. 
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VII.  EXPERIMENTAL  RESULTS 


Experimental  results  for  mechanical  parameters,  as  well  as  resonance  curves  in 
response  to  the  acoustic  and  mechanical  actuation,  are  used  to  refine  the  computer  models 
to  enable  a  more  accurate  and  efficient  design  process. 

A,  OUTPUT  VOLTAGE  MEASUREMENTS 

Two  basic  types  of  measurements  were  performed:  acoustic  and  mechanical. 
Based  on  the  low  output  voltage,  it  was  necessary  to  use  a  Lock-In  Amplifier  for  a 
majority  of  the  measurements.  Further  voltage  and  current  measurements  would  benefit 
from  a  low-voltage  AC  energy  harvesting  conditioning  circuit,  such  as  designed  in  [25]. 

1.  Acoustic  Measurements 

An  SRS-850  DSP  Lock-In  Amplifier  (LIA)  is  used  to  excite  a  Hewlett  Packard 
467 A  Power  Amplifier,  which  drives  an  HL 14-25  Plastic  exponential  horn  (speaker) 
placed  35“  normal  to  the  device;  the  LIA  also  monitors  and  records  the  root-mean-square 
voltage  (Vrms)  output  from  the  device.  In  the  acoustic  set-up  there  is  a  Briiel  &  Kjser 
reference  microphone  with  a  Type  2670  preamplifier  to  measure  sound  pressure  in  the 
immediate  vicinity  of  the  device;  and  a  Polytec  OFV  534  Laser  Unit  which  drives  a 
Mitutoyo  model  K00232104  laser  and  is  controlled  by  a  Polytec  OFV-5000  Vibrometer 
Controller  to  measure  the  maximum  displacement  of  the  device.  Both  the  reference 
microphone  and  the  laser  vibrometer  use  VibSoft  4.7  software  for  data  recording  and 
analysis.  This  set-up  is  shown  in  Figure  20. 
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Figure  20.  Acoustic  test  apparatus. 


Figure  21  shows  the  average  of  Fourier  transforms  of  multiple  acoustic  sweeps 
from  100  Hz  to  1  kHz.  The  data  are  from  two  independent  trials,  each  an  average  of 
5  sweeps.  The  peaks  apparent  at  low  frequencies  are  harmonics  of  60  Hz.  Note  the 
apparent  linear  response  between  output  voltage  and  input  sound  pressure  and  the 
seemingly  independent  displacement  peak  at  522  Hz. 
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Figure  21 .  Acoustic  trials,  average  of  10  sweeps  over  multiple  days, 

100  Hz  to  1  kHz. 


Based  on  the  speaker  manufacturer’s  specifications,  the  reference  microphone 
seems  to  mirror  the  expected  response  curve  from  the  speaker,  seen  in  Figure  22. 


RESPONSE  AND  IMPEDANCE  CURVES  W/  021071  DRIVER  INSIDE  AN 
ANECHOIC  CHAMBER.  1  W  M  m 


•  •  •  lmpe<)«nc«  CofV«. 


Figure  22.  Response  and  impedance  curves  for  HL 14-25  horn,  from  [26]. 
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a. 


Response  to  Varying  Incident  Amplitude 


The  adjustable  power  amplifier  between  the  speaker  and  the  LIA  has  the 
following  gain  settings:  OFF,  xl,  x2,  x5,  and  xlO.  Figure  23  shows  a  linear  response  in 
output  voltage  with  a  ehange  in  incident  speaker  gain.  Also  note  the  increasingly  stable 
phase  angle  as  incident  speaker  gain  gets  higher. 


0  5  10  15  20  25  30  35  40  45  50  55  60  65 


Device  Response  To  Input  Acoustic  Amplitude  at  527  Hz 


Figure  23 .  Linear  output  variation  with  input  gain  adjustment. 
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b.  Modal  Analysis 

The  following  data  was  taken  from  measurements  shown  in  Figure  24.  Analysis 
of  this  data  should  be  sufiieient  to  determine  which  resonance  mode  is  the  dominant 
mode  excited  at  each  frequency. 


Figure  24.  Optical  microscope  image  of  PiezoMUMPs  third  generation 
device,  showing  laser  vibrometer  data  locations. 

Each  leg  is  250  pm  wide. 


Figure  25  shows  the  displacement  per  unit  pressure  vs  actuating  signal  frequency 
obtained  during  multiple  acoustic  sweeps  from  200  Flz  to  900  Flz.  The  data  are  averaged 
from  two  independent  trials,  each  an  average  of  10  sweeps.  The  output  plotted  is 
displacement,  measured  by  the  laser  vibrometer,  divided  by  sound  pressure,  measured  by 
the  reference  microphone. 
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Acoustic  Actuation 
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Figure  25.  Acoustic  trials.  Note  varying  resonance  peaks  at  220  FIz,  525  Hz, 
710  Hz,  and  755  Hz  and  different  relative  sizes  based  on  location 
of  measurement.  The  top  graph  is  at  full  scale  to  show  the  220  Hz 
peak,  the  bottom  graph  is  scaled  to  emphasize  the  higher-frequency 

peaks. 


Figures  26-28  compare  finite-element  models  to  each  of  the  observed  peaks  from 
Figure  25.  The  COMSOL  models  clamp  the  ends  of  the  legs  and  apply  a  normal  uniform 
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pressure  aeross  the  eenter  pad  and  legs.  Note  in  the  measurements  that  the  displaeement 
amplitude  at  the  eenter  eonsistently  measures  substantially  lower  than  the  displaeement 
amplitudes  at  any  other  loeation  aeross  the  frequeney  sweep.  For  the  modal  analysis  the 
effect  of  this  center  measurement  is  considered  qualitatively.  The  one  seeming 
discrepancy  is  in  the  peaks  between  200  Hz  and  250  Hz.  The  likely  cause  is  that  the 
measured  207  Hz  (corresponding  to  COMSOL’s  229  Hz)  peak  has  a  very  low  Q,  so  at 
223  Hz  (COMSOL’s  259  Hz)  a  combination  of  both  is  modes  excited,  thus  the  discrepant 
height  of  the  center  measurement  around  223  Hz. 
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Acoustic  Actuation 


Frequency  [  Hz  ] 


Figure  26.  COMSOL  deformation  shapes  for  modes  at  a)  229  Hz  and  b)  259 
Hz,  eorresponding  to  c)  aeoustie  peaks  observed  at  207  Hz  and  223 
Hz,  respectively.  Note  the  corners  measure  higher  than  the  leg,  and 
the  center  measures  disproportionately  higher  near  207  Hz. 
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Eigenfrequency=579.1  Hz  Surface:  Total  displacement  (a.u.) 
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Figure  27.  a)  COMSOL  deformation  shape  for  mode  at  579  Hz, 

corresponding  to  b)  acoustic  peak  observed  at  526  Hz.  Note  the  leg 
measures  highest,  the  comers  somewhat  in  the  middle,  and  the 
center  barely  deflects  at  all. 
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Acoustic  Actuation 


Frequency  [  Hz  ] 


Figure  28.  COMSOL  deformation  shapes  for  modes  at  a)  710  Hz  and  b)  840 
Hz,  corresponding  to  c)  acoustic  peaks  observed  at  702  Hz  and  751 
Hz,  respectively.  Note  the  leg  measures  highest,  the  comers  do  not 
deflect  at  710  Hz,  corner  1  measures  higher  than  corner  2  at  751 
Hz,  and  the  center  barely  deflects  at  all. 
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c. 


Unusual  Data 


The  223  Hz  peak  appears  to  be  a  linear  combination  of  a  high-Q  223  Hz  peak 
with  a  very  low-Q  205  Hz  peak,  and  is  seen  by  the  laser  vibrometer  as  independent  of 
acoustic  actuation.  This  is  a  much  noisier  signal  than  the  remainder  of  the  Fourier 
transforms,  and  its  origin  remains  unknown.  An  expression  analogous  to  Equation  (8)  can 
be  used  to  find  that  the  approximate  resonance  frequency  of  the  mounting  bracket  is  20 
kHz,  and  laser  vibrometer  measurements  of  the  mounting  bracket  showed  little 
significant  displacement  near  200  Hz. 

The  Helmholtz  effect  was  investigated  of  the  mounting  bracket  that  supports  the 
device,  where  the  device  substrate  acts  as  a  thin  neck  flange.  Viscous  penetration  depth  as 
a  function  of  frequency  is 


S  = 


(13) 


5  3 

where  t]  =  2.0x10'  Pa-s,  and  p  =  1.2  kg/m  are  the  viscosity  and  density  of  air, 
respectively  [27].  At  200  Hz  the  viscous  penetration  depth  is  approximately  150  pm, 
making  the  larger  gap  around  the  outer  rim  of  the  device  (between  the  legs  and  the 
substrate)  the  only  effective  area  of  the  Helmholtz  orifice.  Helmholtz  frequency  for  a 
“thin-neck”  resonator  is 


/o  = 

6  2 

where  c  =  343.6  m/s  is  the  speed  of  sound  in  air,  A  ~  4x10'  m  is  the  total  area  of  the 
effective  orifice,  t  =  10.5  pm  is  the  orifice  thickness  (structure  +  piezo  layers), 
V ~  150x10'  m  is  the  cavity  volume,  and  8  is  the  “end  correction”  and  is  related  to  the 
orifice  area  [27].  Application  of  these  values  to  Equation  (14)  shows  that,  at  200  Hz,  the 
end  correction  8  would  need  to  be  approximately  2  m.  This  is  unlikely. 

Eurther  investigation  of  the  cause  of  this  resonance  is  recommended. 
Additionally,  it  should  be  possible  to  make  a  Helmholtz  cavity  that  stimulates  a  desired 
flexural  mode. 


A 


2^V(t  +  c^)V  ’ 


(14) 
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d.  Lorentzian  Curve-Fitting 

Using  the  above  COMSOL  results  as  a  guide  for  frequeney  and  amplitude  values, 
five  independent  Lorentzian  eurves  were  summed  and  fit  to  the  experimental  data  using  a 
“least  squares”  regression; 


where  Z  is  the  displaeement,  normalized  to  ineident  sound  pressure;  Z  is  a  seale  faetor;^^, 
is  the  eorresponding  peak  frequency;  is  a  measurement  of  the  Full  Width  at  Half 
Maximum  (FWHM);  and  ZBias  is  an  offset  introduced  to  account  for  background  white 
noise. 

The  experimental  data  from  each  point  in  Figure  25  was  averaged;  the  resulting  fit 
of  multiple  Lorentzian  peaks  to  this  averaged  data  is  shown  in  Figure  29  and  summarized 
in  Table  4.  The  quality  factor  is  determined  by  dividing  each  peak  frequency  by  its 
corresponding  FWHM. 
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Average  Displacement 


Frequency  [  Hz  ] 


Leg  Displacement 


Frequency  [  Hz  ] 


Figure  29.  Lorentzian  curves  fit  to  acoustically-actuated  experimental  data 
averaged  over  all  points  (top)  and  only  at  the  leg  (bottom).  The  top 
graph  in  each  set  is  at  full  scale  to  show  the  220  Hz  peak,  the 
bottom  graph  in  each  set  is  scaled  to  emphasize  the  higher- 
frequency  peaks. 
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Table  4.  Summary  of  Lorentzian  Fit-to-Leg  parameters  and  modal  analysis, 
with  relative  displaeement  eolor-eoded  to  mateh  COMSOL 

display. 


Measured 

Freq. 

(Hz) 

Fit  fn 

(Hz) 

Fit 

A/ 

(Hz) 

Fit 

Q 

Center 

Corner 

1 

Corner 

2 

Leg 

Mode 

COMSOL 
Freq.  (Hz) 

207 

210.0 

58.8 

3.6 

High 

High 

High 

Med 

0 

229 

223.0 

223.5 

4.76 

47.0 

Med 

High 

High 

Med 

1 

259 

526.1 

526.9 

7.16 

73.6 

Low 

Med 

Med 

High 

2 

579 

705.2 

707.5 

9.21 

76.8 

Low 

Low 

Low 

High 

3 

710 

752.2 

753.8 

14.6 

51.7 

Low 

Med 

Low 

High 

4 

840 

2.  Mechanical  Measurements 

For  the  meehanieal  vibrational  measurements,  the  same  Loek-In  Amplifier  (LIA) 
and  the  same  energy  harvesting  deviee  were  used.  The  LIA  sends  its  exeitation  signal  to 
an  Aeoustie  Power  Systems  (APS)  Model  114  Dual-Mode  Power  Amplifier,  operated  in 
Current  Mode  for  eonstant  foree  [28],  whieh  drives  an  APS  Model  120S  Perma-Dyne 
Shaker.  On  the  shaker,  in  addition  to  the  device,  is  an  Endevco  Isoshear  Model  770  lA- 
100  Accelerometer.  Both  can  be  either  read  by  the  LIA  or  monitored  by  a  Tektronix  DPO 
2012  Digital  Oscilloscope.  The  device  must  be  sent  through  a  SRS  Model  SR560  Low 
Noise  Preamplifier  in  order  to  produce  a  signal  that  is  clean  enough  to  measure  on  the 
oscilloscope. 
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Figure  30.  Vibrational  test  apparatus. 


A  MotionPro  X3  high-speed  eamera  with  a  Carl  Zeiss  MAKRO-PLANAR  T*, 
2-100  mm,  ZF.2  macro-lens  was  used  to  observe  the  vibrational  mode  at  3910  frames  per 
second  (see  Figure  31).  Software  analysis  of  this  video  will  help  verify  the  dominant 
modes  excited. 


Figure  3 1 .  High  speed  camera  still  image  of  device. 


a.  Output  Voltage 

Figure  32  shows  multiple  resonance  peaks,  the  most  dominant  at  120  Hz  and 
380  Hz,  and  antiresonance  peaks  at  80  Hz  and  360  Hz.  An  analysis  of  the  accelerometer 
output  shows  the  apparent  peak  at  120  Hz  is  artificial;  that  is,  this  harmonic  of  60  Hz  only 
shows  up  due  to  a  dramatic  negative  spike  of  the  accelerometer  output.  Equation  (16) 
shows  the  normalization  method  in  order  to  factor  out  the  variations  of  the  shaker  output. 
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When  necessary,  a  scale  factor  is  used  to  match  magnitudes  to  account  for  minor 
variations  in  output. 


e 


V 

'y^  _  out 


v„. 


—  R  —  f) 

^out  ^  acc 


(16) 


where  Vout  and  Face  are  the  output  voltages  from  the  device  and  the  accelerometer, 
respectively;  similarly  for  phase  angle. 


Mechanical  Actuation 


Frequency  (Hz) 


Figure  32.  Vibrational  sweeps  from  200  Hz  to  900  Hz,  normalized  to 

accelerometer  output. 


b.  Modal  Analysis 

Due  to  the  difficulty  of  including  the  shaker  into  the  acoustic  set-up,  the 
mechanical  measurements  were  taken  in  a  different  room.  This  made  it  impossible  to 
obtain  laser  vibrometer  measurements  to  perform  the  same  sort  of  modal  analysis  as  in 
the  acoustic  measurements;  however,  the  COMSOL-generated  mode  shapes  for  a  given 
acceleration  at  the  leg  bases  are  shown  in  Figure  33.  Note  the  correlation  to  peaks 
observed  at  380  Hz  and  720  Hz,  the  only  modes  where  the  legs  flex  in  opposite 
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directions.  This  suggests  that  only  certain  flexural  modes  will  produce  a  measurable 
potential  difference,  and  that  it  is  opposite  stresses  which  generate  this  potential. 


Figure  33.  COMSOL  mechanically-actuated  mode  shapes  at  a)  366  Hz,  b) 

504  Hz,  c)  661  Hz,  and  d)  800  Hz. 


c.  Lorentzian  Curve-Fitting 

A  more  thorough  analysis  of  the  resonance  peak  at  380  Hz  shows  an  excellent 
Lorentzian  fit,  with  a  high  quality  factor  (0,  as  shown  in  Figure  34.  Two  independent 
measurements  and  a  Lorentzian  fit  are  reflected  in  this  plot.  The  oscilloscope  data  was 
measured  at  discrete  points,  while  maintaining  the  accelerometer  voltage  constant  at 
200  mVnns-  The  LIA  data  was  normalized  to  the  accelerometer  output,  as  above,  and  then 
scaled  using  a  “least  squares”  regression  to  match  the  amplitude  of  the  oscilloscope  data. 
The  Lorentzian  was  then  fit  to  the  LIA  data  using  a  “least-squares”  regression,  from 
377  Hz  -  390  Hz,  ignoring  the  irregularity  at  low  frequencies. 
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Figure  34.  Comparison  of  scaled  LIA  data,  Oscilloscope  data,  and  best-fit 
Lorentzian  at  380  Hz  (top);  showing  phase  shift  on  LIA  (bottom). 


Using  this  fit  data,  Qiarentzian  can  be  determined.  uses  parameters  from 

the  best  fit  Lorentzian  over  the  entire  frequency  range  of  370  Hz  to  390  Hz,  vice  ignoring 
the  irregular  data  at  the  low  frequency,  and  provides  the  most  conservative  estimate  of  Q. 


^Lorentzian 


A/ 


=  106.0 

3.588 


_  /  _  380.861 

^  f^-f_  ^382.314-378.908 


(17) 


^Lorentzian 


380.621 

2.829 


134.5. 


The  380  Hz  peak  is  not  an  artifact  of  the  shaker  or  accelerometer:  first,  since  we 
normalized  our  output  data  for  the  shaker’s  acceleration  (meaning  actuating  force); 
second,  since  the  shaker’s  acceleration  response  levels  out  after  approximately  33  Hz,  as 
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shown  in  Figure  35  and  the  aeeelerometer  has  an  approximately  linear  deviation  (<  2%) 
from  10-1000  Hz,  as  shown  in  Figure  36. 


MAXIMUM  ACCELERATION  LEVELS  FOR  VARIOUS  TEST  LOADS 


Figure  35 .  Acceleration  performance  of  shaker  under  various  test  loads, 

from  [29]. 


TYPICAL  AMPLITUDE  RESPONSE 


Figure  36.  Amplitude  response  for  accelerometer,  from  [30]. 
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3,  Comparison  of  Acoustic  and  Mechanical  Actuation 

Distinctly  different  vibrational  modes  stimulated  by  aeoustie  viee  meehanieal 
exeitation  are  apparent  from  these  measurements.  Due  to  experimental  eonstraints,  some 
possible  eauses  were  not  investigated;  however,  slight  meehanieal  exeitation  in  the 
aeoustie  set-up  allows  for  some  direet  eomparison.  A  Fourier  transform  of  two-minutes 
of  meehanieal  agitation  of  the  aeoustie  experimental  setup  is  shown  in  Figure  37.  Note 
the  223  Hz  peak  is  approximately  5x  larger  in  this  data,  but  does  not  seem  to  show  up 
when  aetuated  by  the  shaker. 
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Mechanical  Actuation 
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Figure  37.  Mechanical  excitation  in  the  acoustic  experimental  set-up.  The  top 
graph  is  at  full  scale  to  show  the  220  Hz  peak,  the  bottom  graph  is 
scaled  to  emphasize  the  higher-frequency  peaks. 


Figure  38  shows  a  comparison  of  this  data  to  the  voltage  data  collected  from 
excitation  using  the  mechanical  (shaker)  set-up.  Note  the  similarity  between 
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corresponding  peaks  and  phase  shifts  at  380  Hz,  600  Hz,  and  approximately  700  Hz.  Also 
note  the  380  Hz  and  600  Hz  peaks  were  not  present  in  the  aeoustie  trials. 


Figure  38.  Comparison  of  meehanieal  aetuation  in  the  meehanieal  and 

aeoustie  experimental  set-ups.  The  measuring  point  in  the  aeoustie 
set-up  is  the  deviee’s  legs  and  has  been  smoothed  for  elarity  of 

eomparison. 


The  COMSOL-generated  mode  shapes  are  repeated  in  Figures  39  and  40.  With 
the  exeeption  of  the  229  Hz  mode,  the  shapes  are  similar  (but  not  identieal)  to  the 
meehanieally-aetuated  mode  shapes.  It  is  likely  the  600  Hz  and  840  Hz  anomalies  in  the 
shaker  data  is  due  to  60  Hz  harmonies,  so  the  only  eorrelating  modes  are  at  380  Hz  and 
aeross  the  two  700  Hz  peaks,  eorresponding  to  opposite  leg  motion.  There  is  too  mueh 
noise  in  the  laser  vibrometer  data  to  perform  mode  shape  analysis  similar  to  the  aeoustie 
trials. 
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Figure  39.  COMSOL  acoustically-actuated  mode  shapes  at  a)  229  Flz,  b)  259 

Hz,  c)  579  Hz,  d)  710  Hz,  and  e)  840  Hz. 


Figure  40.  COMSOL  mechanically-actuated  mode  shapes  at  a)  366  Hz,  b) 

504  Hz,  c)  661  Hz,  and  d)  800  Hz. 
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There  are  only  two  reliable  peaks  shown  on  the  voltage  measurements  of  the 
mechanieally  actuated  device,  381  Hz  and  723  Hz.  A  comparison  of  Figures  38  and  39 
show  the  corresponding  COMSOL  modes  at  366  Hz  and  800  Hz  have  legs  flexing  in  the 
opposite  directions.  Since  the  connecting  pads  are  configured  to  take  the  electric  potential 
from  the  top  of  one  leg  and  the  bottom  of  the  other,  opposite  stresses  are  necessary  to 
generate  a  potential  difference.  Since  the  acoustic  set-up  measures  leg  displacement, 
rather  than  output  voltage,  additional  peaks  are  detected  for  modes  that  do  not  produce  a 
potential  difference. 


Figure  41 .  Lorentzian  curves  fit  to  vibrationally-actuated  experimental  data. 
The  Lorentzian  parameters  for  these  fit  curves  are  listed  in  Table  5. 
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Table  5.  Summary  of  fit  Lorentzian  parameters. 


Measured  Ereq.  (Hz) 

Eit/„  (Hz) 

A/-(Hz) 

Q 

Mode 

COMSOE  Ereq.  (Hz) 

203 

205.5 

8.9 

23.0 

- 

- 

383 

383.6 

7.9 

48.4 

0 

366 

606 

606.2 

10.9 

55.4 

1 

504 

707 

706.2 

7.2 

98.7 

2 

661 

753 

752.9 

8.1 

92.6 

3 

800 

B,  ESTIMATE  OF  DAMPING  PARAMETERS 


Rayleigh  damping  is  eommonly  used  in  finite-element  analysis  and  provides  a 
eonvenient  method  to  eorrelate  the  measured  quality  faetor  to  damping  due  to  the  mass 
and  stiffness  matriees.  The  damping  ratio,  is  related  its  eorresponding  eigenfrequeney, 


4  = 


r  1  ] 

-1- 

do  + 

[2) 

*1  ’ 


(18) 


where  ao  and  ai  are  the  mass  and  stiffness  damping  parameters,  respeetively.  MEMS 
deviees  traditionally  have  low  damping.  In  this  ease  the  damping  ratio  is  related  to  Q 
[31]: 


1 


(19) 


Equating  (18)  to  (19)  and  fitting  the  resultant  equation  to  the  data  in  Table  6, 
initial  trial  damping  parameters  were  determined  to  be  ao  =  22.9  s'\  and  ai  =  2.03x10'^  s. 


57 


Table  6.  Lorentzian  parameters  used  to  fit  initial  Rayleigh  damping 
parameters,  taken  from  the  Lorentzian  curve  fitting 
of  the  acoustic  data. 


Mode 

1 

2 

3 

4 

fn 

223.4 

526.9 

706.9 

755.6 

On 

50.2 

81.2 

97.5 

53.1 

C.  COMPARISON  OF  SIMULATION  AND  MEASUREMENTS 

Finite  Element  Modeling  is  a  useful  technique  for  computer  modeling.  A  smooth, 
continuous  system  can  be  divided  into  a  finite  number  of  discrete  elements,  each  with  its 
own  set  of  differential  equations.  These  elements  share  boundary  conditions,  and  the 
software  numerically  evaluates  the  overall  problem. 

Since  there  was  no  experimental  data  for  basis,  the  COMSOL  finite  element 
models  developed  in  [1]  used  generic,  bulk  material  values  for  Young’s  modulus  and 
Poisson’s  ratio,  as  well  as  selected  damping  parameters  somewhat  arbitrarily.  These 
models  are  further  refined,  and  the  results  described  in  Chapters  VI  and  VII  are  fed  back 
into  the  models  to  achieve  matching  results.  The  COMSOL  models  can  now  be  used  to 
design  the  next-generation  devices. 

1.  Eigenfrequency 

More  careful  and  systematic  selection  of  material  parameters  provides  for  a  more 
realistic  model.  The  most  important  values  that  determine  eigenfrequencies  are  Young’s 
modulus  and  Poisson’s  ratio,  as  they  describe  the  elastic  deformation  of  a  material.  Using 
the  more  accurate  measured  values,  the  COMSOL  generated  eigenfrequencies  are  listed 
in  Table  7. 
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Table  7.  Eigenfrequencies  generated  by  finite  element  modeling. 


Actuation  Method 

Frequency  (f 

iz) 

Center  Pad  (incident  pressure) 

229 

259 

579 

710 

840 

Anchor  (prescribed  acceleration) 

- 

366 

504 

661 

800 

2.  Damping 

Using  the  damping  parameters  calculated  above,  a  frequency  domain  sweep  was 
conducted  in  COMSOL,  shown  in  Figures  42.  Note  the  peaks  are  starting  to  appear 
divided  at  220  Hz  and  710  Hz,  similar  to  the  experimental  data. 


Mechanical  (Anchor  Point)  Actuation 


Frequency  (Hz) 


Figure  42.  COMSOF  frequency  domain  sweep  200  Hz  -  900  Hz  for  center- 
pad  (top)  and  anchor-point  (bottom)  actuation,  in  increments  of 
1  Hz,  using  the  above  damping  parameters. 


The  parameters  entered  into  COMSOF  for  analysis  in  this  thesis  are  compared  to 
the  original  values  and  listed  in  Table  8. 
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Table  8.  Summary  of  COMSOL  input  parameters. 


Young’s 

Modulus  - 

Si  (GPa) 

Poisson’s 

Ratio  -  Si 

Young’s 

Modulus  - 

AIN  (GPa) 

Poisson’s 

Ratio  -  AIN 

ao  (s"') 

ai  (s) 

Original 

170 

0.28 

344 

0.22 

1 

2x10'^ 

Updated 

173.3 

0.28 

384.0 

0.22 

22.9 

2.03x10'^ 

D,  ESTIMATE  OF  POWER  AND  EFFICIENCY 

The  LIA  has  an  input  resistanee  of  50  MO  with  a  25  pF  eapaeitanee,  and  the 
oseilloseope  has  input  resistanee  of  1  MO  with  11.5  pF  eapaeitanee.  Sinee  eapaeitive 
reaetanee  is 


X 


c 


1 

2^’ 


(20) 


the  LIA  values  range  from  63.7  MO  at  100  Hz  to  6.37  MO  at  1  kHz,  and  oseilloseope 
values  range  from  13.8  MO  to  138  MO.  This  provides  an  input  impedanee  range  of  the 
LIA  for  the  frequeneies  analyzed  in  this  thesis  of  50.4  MO  to  81.0  MO,  and  oseilloseope 
impedanee  range  from  13.8  MO  to  138  MO. 

The  maximum  measured  steady  state  voltage  at  resonanee  frequeney  is  6.8  mVp.p 
=  4.8  mVrms  at  380.5  Hz  on  the  oseilloseope  and  220.8  pVrms  at  379.3  Hz  on  the  LIA.  The 
oseilloseope  has  input  impedanee  of  36.4  MO  and  the  LIA  has  input  impedanee  of 
52.7  MO  at  their  respeetive  frequeneies.  This  means  the  deviee  is  providing  P  =  V  /R  = 
0.633  pW  to  the  oseilloseope  and  0.924  fW  to  the  LIA.  Note  the  large  differenee  in 
power  output  for  the  different  monitoring  deviees.  This  superfieially  agrees  with  finite 
element  models  from  [1],  that  there  ean  be  several  orders  of  magnitude  differenee  in 
power  output  depending  on  load  impedanee. 

Of  note,  these  values  are  power  delivered  to  equipment  that  is  designed  not  to 
draw  mueh  power.  Matehing  impedanees  using  a  eireuit  sueh  as  designed  in  [25]  will 
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enable  an  inerease  in  power  produced  per  device  of  several  additional  orders  of 
magnitude.  A  summary  of  this  estimation  is  below  in  Table  9. 


Table  9.  Power  delivered  by  one  device. 


Voltage 

(rnVrms) 

Frequency 

(Hz) 

Impedance 

(MO) 

Power 

(iW) 

Oscilloscope 

4.8 

380.5 

36.4 

633.0 

Lock-In  Amplifier 

0.221 

379.3 

52.7 

0.924 
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VIII.  FUTURE  WORK 


Most  of  the  Future  Work  reeommendations  from  [1]  and  [25]  remain  valid.  These 
inelude  refinement  of  the  mierofabrieation  proeess  and  analysis  of  eombinations  of 
eonditioning  eireuitry  for  multiple  deviees.  Further  reeommendations  are  listed  below. 

A,  MODELING 

With  the  data  obtained  in  ehapters  VI  and  VII,  more  detailed  and  aeeurate 
eomputer  models  ean  be  obtained.  These  models  ean  then  be  used  to  design  better,  more 
optimal  deviees. 

1.  Exploring  the  Parameter  Space 

There  are  several  parameters  that  would  benefit  from  further  analysis,  and  more 
aeeurate  matehing  of  finite-element  models  to  experimental  values.  Eigenfrequeney  and 
modal  matehing,  as  well  as  quality  faetor  matehing,  are  the  most  important  of  them. 

Comparison  of  the  high-speed  video  deseribed  in  Chapter  VII,  and  of  additional 
videos  at  other  frequeneies,  to  COMSOL-generated  mode  shapes  will  enable 
eonfirmation  of  the  dominant  modes  exeited  by  meehanieal  stimulation.  Further 
investigation  of  laser  vibrometer  results  at  various  loeations  on  the  deviee  and  their 
eomparison  to  COMSOL  mode  shapes  will  also  enable  eonfirmation  of  the  dominant 
modes  exeited  by  aeoustie  stimulation. 

Continued  analysis  of  material  (Young’s  modulus  and  Poisson’s  ratio,  for 
example)  and  damping  parameters  (mass  and  struetural  damping)  is  neeessary  to  further 
mateh  finite-element  models  with  experimental  results.  More  realistie  and  reliable 
eomputer  models  ean  then  be  used  to  test  future  generations  without  requiring  as  mueh 
experimental  validation.  This  will  signifieantly  reduee  the  timeline  and  the  eost  of  design. 

2,  Next  Generation  Devices 

Other  than  the  obvious  adjustments  of  leg  and  eenter  pad  dimensions  and 
geometry,  the  next  generation  deviees  eould  be  improved  by  determining  the  points  of 
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highest  potential  differenee  and  laying  the  eleetrode  pads  aeross  the  equipotential 
surfaees.  One  option  is  to  eover  the  top  of  the  entire  piezo  layer  with  a  single  A1  pad 
metal  layer  and  eonneet  the  bottom  of  eaeh  piezo  leg  together  with  another  pad  metal 
layer  to  obtain  the  highest  possible  eurrent  for  the  “top-to-bottom”  potential  differenee. 

B,  FABRICATION 

The  ability  to  deposit  alternative  strueture  layers  and  piezo  layers  at  NFS  will 
greatly  reduee  the  eost  of  eaeh  MEMS  deviee. 

1.  NFS  Fabrication  with  SOI  and  SU-8 

With  the  Boseh  deep  reaetive  ion  eteh  instrument  fully  funetional,  it  is  now 
possible  to  finish  mierofabrioation  of  the  designs  intended  for  NFS  fabrieation.  Final 
refinement  of  the  fabrieation  proeess  and  analysis  of  SU-8  as  a  strueture  layer  will  be 
neeessary,  as  listed  below. 

a.  Fabrication 

•  Determine  the  Boseh  etehing  parameters  for  baekside  trenehing  to  free  a 
suspended  deviee. 

b.  SU-8  Analysis 

•  Determine  experimental  data  for  deposited  SU-8,  sueh  as  Young’s 
modulus,  hardness,  and  Foisson’s  ratio  to  update  eomputer  simulations. 

•  Determine  eigenfrequeneies  of  eaeh  multileg  strueture  and  eompare  to 
eomputer  simulations. 

•  Determine  resilienee  of  SU-8  struetures,  both  long-term  eyelie  fatigue  and 
one-time  failure. 

2.  NFS  Deposition  of  AIN 

•  Continue  analysis  of  reaetive  sputtering  AIN  at  NFS. 

•  Investigate  other  options  for  deposition  of  AIN,  sueh  as  reaetive  ion 
impregnation  with  a  plasma  eteher. 
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C.  TESTING 

Initial  testing  of  a  single  deviee  was  eondueted  in  this  thesis.  More  extensive 
testing  is  neeessary  to  refine  finite-element  modeling  parameters  to  input  into  COMSOL. 

1,  Conditioning  Circuit 

Sinee  the  power  output  from  eaeh  deviee  is  a  funetion  of  the  load,  determining 
output  impedanee  of  eaeh  deviee  is  neeessary  to  design  the  optimum  eonditioning  eireuit. 
Very  little  modifieation  to  the  eireuitry  designed  in  [25]  will  be  required  for  this  purpose. 

2,  Multiple  Device  Arrays 

Since  different  resonance  modes  will  generate  outputs  with  different  phases,  it  is 
possible  that  the  optimal  arrangement  of  multiple  devices  will  differ  based  on  the 
dominant  mode  excited.  An  in-depth  analysis  of  output  waveforms  in  conjunction  with 
utilization  of  the  conditioning  circuit  above  will  provide  the  first  steps  toward  developing 
an  array  of  MEMS  vibration  and  acoustic  energy  harvesters  that  can  provide  usable 
power. 


3.  Power  and  Efficiency  Analysis 

Updated  COMSOL  data  will  enable  better  estimates  of  the  output  power  and 
efficiency  of  individual  devices  and  of  arrays  of  devices,  which  can  be  compared  to 
measured  values.  Following  this  determination,  a  Ragone  plot  showing  power  density  vs. 
energy  density  will  allow  a  convenient  comparison  to  other  energy  harvesting  devices. 

4.  Failure  Analysis 

The  device  analyzed  in  Chapter  VII  has  been  subjected  to  several  hours  of  use 
with  no  discemable  change  in  behavior.  MEMS  devices  are  characteristically  robust; 
however,  a  deeper  understanding  of  the  effective  lifetime  of  any  device  can  be 
extraordinarily  helpful. 

A  small  selection  from  each  design  must  be  subjected  to  extreme  conditions  for 
failure  analysis,  and  the  rest  must  be  tracked  over  their  lifetime  to  monitor  fatigue.  Useful 
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information  to  track  would  be  the  effeet  of  numerous  eyeles,  normalized  to  aeeeleration 
or  pressure,  on  eigenfrequeney  values  and  peak  displaeement  and  output  voltage. 

5,  Investigating  Unusual  Findings 

During  Young’s  modulus  determination  unusual  data  was  diseovered  on  the 
surfaee  (50  nm-100  nm)  of  the  AIN  layer.  Further  analysis  of  the  surfaee  of  the  AIN 
layer  will  help  build  even  more  aeeurate  eomputer  models. 

A  resonanee  peak  between  200  Hz  and  250  Hz  was  deteeted,  even  without 
ineident  sound,  during  the  aeoustie  trials.  COMSOL  models  prediet  the  two  lowest 
resonanee  frequeneies  to  be  near  this  value.  Investigation  into  the  eause  of  this  peak  will 
lead  to  a  better  understanding  of  the  MEMS  device  and  its  effeets. 

D,  MODAL  ANALYSIS 

The  only  resonanee  frequeney  that  seems  to  overlap  between  experiments  using 
aeoustie  and  meehanieal  aetuation  is  between  700  Hz  and  750  Hz.  The  COMSOL- 
generated  mode  shapes  at  840  Hz  (aeoustie)  and  800  Hz  (meehanieal)  are  very  similar. 

1.  Center-Pad  versus  Anchor-Point  Actuation 

The  resonanee  frequeneies  are  different,  both  in  value  and  in  mode  shape,  for 
different  methods  of  aetuation.  Analysis  of  the  differing  effeetive  masses  for  eaeh  method 
of  aetuation  will  enable  understanding  of  this  phenomenon. 

2.  Helmholtz  Resonance 

Voltage  is  generated  when  the  deviee  legs  flex  in  opposite  direetions.  Aeeording 
to  COMSOL  this  oeeurs  at  259  Hz,  579  Hz,  and  840  Hz  for  eenter-pad  (aeoustie) 
aetuation;  and  at  366  Hz  and  800  Hz  for  anehor-point  (meehanieal)  aetuation.  Experiment 
shows  a  possible  voltage  peak  near  700  Hz  for  aeoustie  aetuation.  It  should  be  possible  to 
design  a  eavity  that  stimulates  this  resonanee  frequeney  to  enable  harvesting  of  ambient 
aeoustie  energy. 
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